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Solutions of the Einstein equations which are periodic and have standing gravitational waves, in 
the weak-field zone, are valuable approximations to more physically realistic solutions with outgoing 
waves. A variational principle for the periodic solutions is found which has the power to provide, for 
binary systems with weak gravitational radiation, an accurate estimate of the relationship between 
the mass and angular momentum of the system, the masses and angular momenta of the components, 
the rotational frequency of the frame of reference in which the system is periodic, the frequency 
of the periodicity of the system, and the amplitude and phase of each multipole component of 
gravitational radiation. Examination of the boundary terms of the variational principle leads to 
definitions of the effective mass and effective angular momentum of a periodic geometry which 
capture the concepts of mass and angular momentum of the source alone with no contribution from 
the gravitational radiation. These effective quantities are surface integrals in the weak-field zone 
which are independent of the surface over which they are evaluated, through second order in the 
deviations of the metric from flat space. The variational principle provides a powerful method to 
examine the evolution of, say, a binary black hole system from the time when the holes are far apart, 
through the stage of slow evolution caused by gravitational radiation reaction, up until the moment 
when the radiation reaction timescale is comparable to the dynamical timescale. 
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^ ■ I. INTRODUCTION 

o ' 

The difficulties associated with the study of solutions of the Einstein equations are well known. However, a number 
of specialized techniques have been developed which allow analyses in rather narrow circumstances. For example, 
the post-Newtonian approximation yields relativistic corrections for systems which are restricted to slow speeds and 
■ — , , weak gravitational fields. Or, symmetric geometries with two or three Killing vectors allow for study of Kerr or 
^ ' Schwarzschild black holes and both rotating and non-rotating neutron star models. A particularly fruitful technique 
I \ involves the perturbation analysis of analytically known solutions — perturbations of flat space comprise linearized 
£h ■ gravity H ; the perturbations of spherically symmetric geometries are used to study the emission of radiation from 

test particles orbiting a black hole and the quasi-normal oscillations of black holes || and neutron stars || . 
^ . This paper presents a new restriction of the Einstein equations which is sufficiently limiting that analysis can 
proceed and yet sufficiently general to encompass a wide variety of interesting applications. Interest focuses on those 
solutions of the Einstein equations which allow a coordinate system in which the geometry is periodic in time inside 
a bounded region of space-time. Such a geometry might have strong fields, gravitational waves and high speeds; and 
it might involve black holes or neutron stars. 

In quantum mechanics the physically realistic solutions of the Schrodinger equation with, say, an incoming wave 
packet being scattered into an outgoing wave packet are constructed from a linear combination of the periodic solutions 
which contain standing waves. Thus, the scattering problem is reduced to the analysis of periodic solutions. 

In general relativity the nonlinearity of the Einstein equations intrudes on using this same idea of constructing 
physically interesting solutions from the periodic ones. But with care and some limitations, the basic process still 
provides an avenue toward otherwise unapproachable physical systems. 

Periodic solutions of the Einstein equations with radiation are not asymptotically flat || — this is not surprising: 
in the weak-field zone, where the linearized Einstein equations give an approximate description of the gravitational 
field, the mass density of standing waves falls off as r~ 2 . Thus, the contribution of radiation to the mass content 
inside a radius r grows linearly with r; and the linearized Einstein equations cannot give an accurate description of 
the gravitational field out to either spatial or null infinity. But this is not a critical limitation to periodic geometries. 

In this paper much attention is focused on the "weak-field zone" — a region of space-time, surrounding a fully 
relativistic source, throughout which the linearized Einstein equations give an accurate description of the gravitational 
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field. A further requirement which we impose upon our use of the weak-field zone is that its total energy content in 
gravitational waves be much less than the mass of the source. Thus the weak-field zone is restricted to a region where 
the nonlinear effects of the radiation are small. For modeling realistic binary systems the boundary of our region 
of interest is always within the weak-field zone. Thus, for periodic geometries we consider only a bounded region of 
space, and describe and impose boundary conditions well away from spatial or null infinity — which might not even 
exist for some continuations past the boundary of a periodic geometry! 

For a periodic example, consider the binary pulsar — the neutron stars themselves are relativistic, collapsed objects; 
and the system has an orbital period of a bit less than eight hours and a total mass of about 3M Q . A region bounded 
at a distance of one light month contains approximately 200 wavelengths of quadrupole radiation from the source; 
m/r ~ 10~ n ; the approximate energy content in the gravitational waves out to this distance is 1O _22 M0; and the 
approximate binding energy of the system is Gm 2 / R ~ 7 x 1O~ 7 M0. The true, physical system has outgoing radiation 
near the boundary, and the orbit decays as a consequence. But a similar bounded, periodic, exact solution to the 
Einstein equations could be constructed from this physically realistic one by sending gravitational radiation inward 
from the boundary with amplitude and phase chosen to keep the system from evolving. The nature of the periodic 
geometry outside the boundary need not be considered. However, if the periodic geometry were extended outward 
to a distance of about 10 22 light years, then the energy of the radiation would dominate that of the source; and, 
even at linearized order and with ignorance of the time evolution, the periodic geometry would no longer resemble 
the physically realistic one. For this system, only at such a large distance does the lack of asymptotic flatness of the 
periodic geometry intrude upon the analysis. 

A second example is provided by a test particle of small mass [x close to the innermost stable circular orbit about 
a black hole of mass M with /i <C M. A perturbation analysis, shows that the test particle nearly moves along a 
geodesic; the secular corrections are from radiation reaction effects at order (/i/M) 2 . The physically reasonable solution 
to this system involves only outgoing radiation. But a corresponding exact periodic solution can be constructed by 
demanding standing wave boundary conditions at a large but finite distance from the black hole. The energy density 
of the gravitational waves at a distance r is proportional to (n/Mr) 2 . As long as the radius r of the boundary is large 
enough that r 3> M but small enough that r{n/M) 2 < M then the energy content of the radiation will be dwarfed 
by the mass of the black hole, and the boundary will still be in the weak-field zone. 

These periodic solutions of the Einstein equations are interesting, involve strong gravitational fields, high speeds, 
are not asymptotically flat and cannot be described by the linearized Einstein equations except in the weak-field zone. 
These are typical of the solutions of interest in this paper. 

The relationship between periodic geometries and physically realistic ones with outgoing waves was carefully con- 
sidered in Paper I || . There it was shown that a specific linear combination of similar periodic geometries, of differing 
frequencies with each geometry containing only standing waves in the weak-field zone, was an approximation to an 
exact solution of the Einstein equations with outgoing waves. Similar methods were first used in the context of general 
relativity by Thorne H who analyzed perturbations of neutron star models. 

The error in this approximation was also carefully analyzed in Paper I. When the effect of radiation reaction is 
weak, the linear combination is sharply peaked at a resonant frequency. And some physical quantity, describing an 
aspect of the linear combination of metrics, differs from the corresponding quantity for the exact outgoing-wave metric 
by an amount comparable to its change in one cycle due to the effects of radiation reaction. Thus, a resonant periodic 
solution alone accurately models an outgoing-wave solution as long as the radiation reaction timescale is much longer 
than the dynamical timescale. And the error in the linear combination is inversely proportional to the ratio of these 
two timescales; this is also proportional to the ratio of the frequency width of the linear combination to the resonant 
frequency. 

Relativistic binary systems are, perhaps, the most interesting of the systems which can be approximated in this 
way. In particular the gravitational wave luminosity and the location of the innermost stable orbit in a relativistic, 
binary system can be studied. Other phenomena of interest which are approachable via periodic geometries include 
axisymmetric, rapidly rotating neutron star models and quasi-normal oscillations of black holes or neutron stars. 
However, even in circumstances when radiation reaction forces are strong and evolution is normally rapid, the periodic 
solutions are still interesting in their own right — after all they are solutions of the Einstein equations with generally 
strong gravitational fields and gravitational radiation. 

Thus far the two most useful methods of studying relativistic binaries are the post-Newtonian approximation and 
the test particle approximation. In both of these approaches the effects of radiation reaction are ignored at the first 
order of a small quantity, and the consequent orbits studied are conservative and periodic. The periodic assumption 
includes both of these approaches as special cases and is, therefore, of more general validity. 

Paper I contains a variational principle restricted to time independent geometries; in this paper the variational 
principle is generalized to include periodic geometries. The generalization reveals the role of gravitational radiation 
more clearly and substantially simplifies the treatment of the radiative boundary terms in the weak-field zone. Also, 
in Paper I the boundary terms had to be evaluated in the local wave zone; in this paper it is necessary to go out only 
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to the weak-field zone where the geometry is accurately approximated by the linearized Einstein equations — for the 
Sun-Earth system this is the difference between going out to Alpha Centauri for the local wave zone and out to the 
orbit of Jupiter for the weak field zone. And given an approximate periodic geometry, which differs from an exact 
solution to the Einstein equations by order (5, the variational principle has the power to yield accurate estimates, 
with an error of order S 2 , of all the interesting quantities which describe the geometry, except the actual metric itself. 
These accurately estimated quantities include the effective mass and angular momentum of the total system as well 
as its individual components, the frequency of periodicity, the angular frequency of rotation, and the amplitudes and 
phases of the gravitational radiation in a multipole decomposition. 

The quantity being extremized in the variational principle is called the effective mass of the system. It is defined as 
a surface integral in the weak-field zone, and has the interesting property of being independent of the actual surface 
chosen for the integral, through terms of second order in the deviation of the metric from flat space as long as the 
surface is in the weak-field zone. In particular, the effective mass does not include the mass content of the standing 
gravitational waves, which diverges and keeps the geometry from being asymptotically flat. And a second surface 
integral, involved in the variational principle, defines the effective angular momentum in a similar manner which also 
does not include a divergent contribution from the standing gravitational waves. 

Section |n| first reviews the 3+1 formalism of general relativity and, then, reviews and modestly modifies Thome's 
unique notation for symmetric trace free tensors, outlines his general solution of the linearized Einstein equations, 
and transforms to a gauge more suitable for present purposes. Previous familiarity with Thome's notation greatly 
facilitates the reading of this paper. Section III gives a variational principle for a generic, periodic solution to the 
Einstein equations. The actual quantity being extremized is closely related to the quasi-local energy as developed by 
Brown and York Q , and the appropriate boundary conditions for the variational principle include the specification of 



the three metric on the sides of the bounded region of space-time under consideration. In §IV the variational principle 
is modified for the special case that the boundary is in the weak-field zone, and the boundary integrals are rewritten 
in terms of the amplitude and phase of the gravitational waves of different multipoles. The variational principle 
leads to formal definitions of the effective mass and effective angular momentum of a periodic geometry. Some of the 
complicated analyses are described in the Appendix. 



II. BACKGROUND AND NOTATION 



A. The Initial Value Formalism and Dynamics of General Relativity 

A four dimensional space-time with a metric g ab may be foliated into spacelike hypersurfaces of constant t, with a 
metric 7 ab , §] 

g ab dx a dx b = -N 2 dt 2 + -y ab (dx a + N a dt){dx b + N b dt). (1) 

The quantity N is the lapse function, and N a is the shift vector. The three dimensional metric has a derivative 
operator D a , and Ricci tensor R ab . 

The extrinsic curvature, K ab , of the hypersurface is defined from 

Q ab = -2NK ab + 2D (a N b) = £ tlab , (2) 

where Ct is the Lie derivative with respect to the time translation vector, t a d/dx a = d/dt, which points in the 
direction of increasing t with all spatial coordinates held fixed. In the Hamiltonian formulation of general relativity 
, the momentum conjugate to 7 a (, is 7r ab /167r where 

ir ab = -7 1 / 2 (K ab - ~/ ab K c c ) , (3) 

and 7 is the determinant of j ab . 

The constraint equations on a given hypersurface are restrictions on -f ab and Tr ab from the Einstein equations. These 
are the Hamiltonian constraint, 

N = R + 7~ 1 (^a < V - Tr ab ir ab ) = 16np, (4) 

and the momentum constraint, 

N a = D^h 1 ' 2 ) = ~8Trj a , (5) 
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where p is the energy density and j a is the momentum density of the stress energy of matter, 

p = N 2 T U , (6) 
f = Nl ab T b \ (7) 

and the spatial part of the stress-energy tensor, 

Sab = la C lb d T cd , (8) 

is used below. 

The dynamical part of the Einstein equations gives 



Tjab — 7\r 1/2 f jjab ^ „ab r> \ i Ar~, — l/2„,ab ( 'cd^ ^ c d 

P = -iV7 ' IK - -7 K J + -iV7 '7 I 7r 7r C(i - -7r c 7r fl 

-2A^ 7 - 1 / 2 (V c vr c b - ^ C c7r ab ^) + 7 1/2 (D a D b N - y ib D c D c N) 
W /2 D C (7" 1/2 7r ab 7V c ) - n ac D c N b - Tr bc D c N a + 8nN-f 1/2 (S ab - P1 ab ) 



C t TT at >. (9) 



B. The Linearized Analysis in the Weak-Field Zone 



The space-times of interest in this paper are periodic and may contain gravitational standing waves far from the 
source. And periodic, standing waves cannot extend to spatial infinity in an asymptotically flat geometry However, 
if the amplitude of the waves is small enough then a weak-field zone around the source exists wherein the geometry is 
well approximated by the linearized Einstein equations. It is in this weak-field zone that we analyze the gravitational 
waves and define quantities similar to the mass and angular momentum of the system. 



1. Notation 



For the geometry in the weak-field zone, a slight modification is made of Thome's notation for coordinates and 
symmetric trace- free tensors. The lapse function at zeroth order in the deviation from flat space is a constant Nq, 
not necessarily unity; thus the Minkowskii time coordinate is Not, and the Cartesian coordinates (x,y,z) are the 
traditional flat space coordinates, with (r, 9, 4>) being the corresponding spherical coordinates. Indices (i, j, k, p, q, r, s) 
run over (x, y, z) and denote the Cartesian components of a spatial tensor with a flat metric. Summation is implied 
when these indices are repeated, and such indices are sometimes both lowered without ambiguity. The flat space 
derivative operator is Vj. A comma denotes a partial derivative with respect to a Cartesian coordinate. And the 
quantity = V^r is the unit outward radial vector. 

Tensors with large numbers (say I) of indices are common. Thus in the convenient abbreviation 

= Thh-i, (10) 

the indices represent tensor components in a Cartesian coordinate system, just like k,p, q, r, s) do. Furthermore, 
often such tensors are both symmetric and trace-free (STF) on all pairs of indices and also have Cartesian components 
which are functions only of t and r and independent of 9 and <j); the Cartesian components of such STF tensors are 
written as capital script letters. For example, 

Tjihh—h = (12) 

and 

dT Jl /d6 = 0, dTjJd(t> = Q. (13) 
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Also the abbreviation 

Rk, =r kl r k2 ...r kl (14) 

is useful to represent the outer product of many unit radial vectors. 

A convenient decomposition of the STF Z-tcnsors is in terms of a basis set of 21 + 1 constant STF tensors y 1 ™, with 
—l<m<l, defined by Thornc (his Eq. (2.12)). Two useful properties of the y l ™ are 

Y lm =y l j™R Jn (15) 
where Y lrn is the usual spherical harmonic function, and their orthogonality 



llmylm'* = &l + 1)U 

With this basis set of STF tensors the decomposition of Tj l (t, r) yields 

i 

T Jl (t,r)= TUt,r)yf. (17) 

m— — I 

Also if Tj t (t,r)/r is a solution to the flat space wave equation then it may be separated into outgoing and ingoing 
parts, 

T Jt (t, r) = T%*(N t -r)+ T™{N t + r). (18) 

And the decompositions, 



T^(N t - r) = £ T^(N t - r)yf, (19) 

m— — l 
I 

Tj«(N t + r) = ]T T Z(Not + r)yf, (20) 



ra— — l 



are natural. Furthermore, if and are each composed of a sum of periodic pieces of different frequencies u mn , 
then complex amplitudes and phases 9f^ are defined from 

T%?(N t -r)=J2 TL exp[-iC + + ^mn{N t - r)], (21) 

n 

and 

TZ (Not + r)=Y^ exp^L" - + ^ mn (N t + r)} . (22) 

n 

If the waves are standing, then the 9j™ are real so that the outgoing and ingoing magnitudes are equal. The above 
equations lead to 

T lm (t, r) = ]T ^Tm cos(C" - \l* + ^)^* (23) 

n 

and 

1 1 

T * (*. r ) = E E 2T "^ cos (C - 2^ + ^nOe--^*. (24) 

n m— — l 

Thorne defines yj 71 as a special STF tensor which depends upon a particular orientation of the Cartesian axes. 
Thus one frame of reference rotated with respect to a second has a different corresponding y™. For a frame of 
reference rotating with a coordinate angular velocity about the z-axis with respect to an inertial frame 
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y Jl rot _ 6 ^Jl inortial' \ Z0 > 

from Thome's Eq. (2.12) if the axes are aligned at t = 0. And as viewed from the rotating frame 

1 1 
T * (*. ^EE 2T /™^r rot cos(C " 2^ + <v)e i(aw +™n/w )Wo*. (26 ) 
n m=— i 

The final specialization of interest is that 7j, be periodic with fundamental frequency wo as viewed from the rotating 
frame of reference. Eq. ( |26| ) then requires that 

Wmn = nuo - mCl/N 0) (27) 

for n an integer. 

If 7j ( is constant in time in the inertial frame of reference then uj mn must vanish for all m and n for which T™ m is 
not zero. If, in addition, is nonzero and Q/uiqNo is not rational then the constant = unless m = n = 0, i.e. 
T/ ; is constant in time and axisymmetric; or if f2 is zero then TJ™ = for n 0. Also, if Q/ljoNq is an integer then 
T/ ; is also periodic when viewed from an inertial frame of reference. 

From Thome's definition of y 1 ™ it follows that 

yhn = ^mylj-m^ (2g) 

and consequently if Tj l is real then 



q T,n 
7 Z,m 



nT,-n 



(30) 



and 



T? m = (-1) 



l-\-jnj-i — n 



(31) 



,2. General Linearized Solution 



Thorne Q gives a general solution to the linearized Einstein equations in terms of the mass and current moments 
evaluated at a retarded time, Ik ( {Not — r) and Sk , (Afci — r), in one specialization of the Lorentz gauge. It is convenient 
to separate the stationary, time independent parts of the metric from the (often radiative) time dependent parts. The 
superscript refers to the constant moments so that 1° and S® are the constant mass monopole and current dipolc 
moment of the geometry — the mass and angular momentum of the geometry if there is no gravitational radiation. 
Useful definitions are 

2X° ~ 2(2/ -1)!! ^, 



i=i 



and 



_ -2p*e lpq S°r q " 4Z(2Z-1)!! 



2=2 



The boldface €j pq is the Levi-Civita tensor. 



Thome's metric is just a perturbation away from flat space-time, and consequences of his Eq. (8.13), through first 
order in the perturbation, are 



{-. DO , -d 

i=2 



(34) 



and 



G 



N* = N S> + N f ji J2 n + 1)[ [r-'e^S^ (N t - er)] 



1 = 2 



oc . , w 

-N r £ (Not - er)], Kl _ t 



1 = 2 



(35) 



for the lapse and shift, a dot represents a derivative with respect to the function argument (Not — er) and e is +1 for 
outgoing radiation and —1 for ingoing. The three- metric, through first order, has the stationary part separated from 
the time dependent part by 

7ij = fa(l + I) + hij (36) 
where is the fiat three metric in Cartesian components, and the time dependent terms are all contained in 



°° 2f-eV 



1=2 
oo 

E 

1=2 
oo 

E 

1=2 



4(-e) 



i r 



r 1 i ijKl _ 2 (N t - er) 



,K, 



8l(-e 
T + T)! 



r 1 e. pq (iS J ) pKl _ 2 (N t - er) 



(37) 



A different version of the Lorentz gauge is preferred here wherein all of the time dependence is removed from the 
lapse and shift vector. The gauge change generated by (£*,&), where 



d&fdt = -A 2 E V-^K (Not - er)] , Kl 



1=2 



dij/dt = -Vj-ft - (Nj - N Sj), 



(38) 



accomplishes this. The ultimate result is that the general solution to the linearized Einstein equations, in the preferred 
gauge, has a lapse and shift vector 



and 

And jij is still given by Eq. ( |36|) but with 

,2 



N = N (l — 1/2), 
N 3 = N S j . 



(39) 
(40) 



flin 



E(- £ ) Z 7j {Mrt^U +2[r- 1 xg f 2 i ) _J,^_ 2 +[r- 1 I { K ; 2 \ l3Kl -4[r- 1 X Ki _ l(l ], j)K; _ 1 } 



1=2 



\ - (~e) l+1 8l r _ x (+ i) i 1 

/ , (l + 1)1 V r e Pl( i0 3)pK l _ 2 \^ K i-2 L~ ^piCi-j^iJ'i)^!-! J • 



(41) 



where the parenthesized superscript, e.g. Spjc? 1 ; denotes differentiation, with respect to the implied argument 
(Not — er) the appropriate number of times, positive or negative. Also 



1 dhij 



In this gauge a number of useful identities hold: 



(42) 



(43) 



r 3 K» = o, 



(44) 
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= 0, (45) 

nS l = 0, (46) 

V,S l = 0, (47) 

V fc V fe 7 = 0, (48) 

W k W k S l = 0, (49) 



and 

1 d 2 h. 

In the local wave zone the leading 1/r contributions to hij from the time dependent moments give 



V,V fc ^-^^f = 0. (50) 



1=2 

°° 81 1 
where troy is the two dimensional metric of a constant-^ two sphere. Also in the wave zone r l hij = (D(r~ 2 ) 



(51) 



3. The Standing Wave Solutions in the Rotating Frame of Reference 

The amplitudes If m and phases 6{™ of the gravitational radiation are defined in an inertial frame of reference in a 
manner similar to Eqs. and (§|): 

{-e) l l Kl (N t - er) = £ exp[-ie(0£ - -Jtt) + i^w^o* - er)] (52) 

and similarly for Sk, ■ 

For the special case that the geometr y is pe riodic in a rotating frame of reference, and contains no traveling waves 



(only standing waves) the notation of § 1IB 1 allows the linearized geometry in the rotating frame of reference to be 
written in terms of the stationary moments, lf m and S? m , and the periodic, radiative moments i|L and S™ m as 



and 



where 



and 



JV = JV (l-7/2), (53) 

= N„S j + n& (54) 

&d/dx l = d/dcj), (55) 
2J° V 2(2/-1)!! . 

1 - — + r l+ll\ h^Kx R K, , (56) 



_ 2p l e ipq r q S° \ - 4/(2? - 1)!! ^ im 



8 



Also 



l 7 rn,n 

-2u 2 mn y l ^ Kl _ 2 [r- 1 cos(C ~ ^ + u> mn r)] >Kl _ 2 

-^ly^ir' 1 cos(0g£ - -Ztt + uw)],^, -4^_ l(i [r" 1 cos(C - \^ + fW)]>i)jr,_i 
£ 77^^™°*°' {"mne^y^Jr- 1 sin(C - ^ + u mn r)], qKl _ 2 

+u m 1 n y l p K l - 1 £ Pq ( 1 [r' 1 sin(C - ^Ztt + cw)].,-)**,^ | . (58) 



In the wave zone this becomes 

I 



/, 8(w m „) 7-n „„„/"/)nJ , ,, ..MiUm 73 /„P „9 ^„.J>9„. \AnuoNot 

hi i 71 I im cos ^im+ UJ mnr)y Kl _ 2Pq KK l _ 2 {cr 0l cr 0:j - -a a 0lj )e 



+ r(7 + 1)! ' m ^ ' m " fw ™" r ) ers pVi-2 r R K l _ 2 (cr 0i a 0j - -a Q <j Qlj )e . (59) 

These equations give the general periodic solution to the linearized Einstein equations, in the gauge described above, 
with stationary moments, lf m and Sf m , and constant wave amplitudes and phases, {I^ m , #f") and (S™ m ,9f£). And for 
hij to represent real, standing waves it is necessary that 

iL = (-i) l+m ir™, (so) 



9 In 
Im 



nI,—Tl 

'L—m'> 



(61) 



and similarly for S™ m and 6f£. 



Finally, it is useful in § IV and in the Appendix to separate out the part of hij in Eq. (p8f) which includes contributions 



from both If m and and is summed over I, with both m and n held fixed. Thus 

h - = \ ~* hV %n 



(62) 



and 



Vntimn 2 jmn r, 

fcV fty - ^> m n h ij - U - 



(63) 



III. A VARIATIONAL PRINCIPLE FOR PERIODIC SOLUTIONS OF THE EINSTEIN EQUATIONS 

A variational principle comes from the traditional Hamiltonian formalism of general relativity and is very closely 
related to the Hamiltonian and the concept of quasi-local energy as developed by Brown and York [Q, their Eq. (4.13). 
The starting point is the definition 



167T#l 



N 



+ 2N a D b (n b a / 1 1 / 2 )}^d 3 x 
j> 2N b 1 ~ 1/2 n b a r a ^d 2 x- j> 2Na ab D a r b y/ad 2 x. 

J rn J rn 



(64) 



This volume integral is over a spacelike hypersurface of finite extent which is bounded by a two-surface defined by a 
scalar field, r, which is constant on the boundary, r = tb- Here in §11] the vector r a is the outward pointing unit 
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normal to the bounding two-surface, a ab is both the metric of the two-surface and the projection operator onto the 
two-surface, 

Cat = lab - r a r b , (65) 

and a is its determinant; and there is no restriction on the location of the boundary of the finite region in which we 
are interested — in particular it is necessary neither that the boundary extend out to the weak-field zone nor that the 
scalar field r which defines the boundary be related to a flat space radial coordinate. 

An arbitrary, infinitesimal variation of N,N a ,-f ab and ir ab , which holds the location of the boundary fixed, results 
in an infinitesimal change in Hi , 

16nSH 1 = - J (SNATyfy + 28N a N a ^ + H,bV ab - 8ir ab G ab ) d 3 x 

+ <f (r a N a ~{- l ' 2 ir bc 8~f bc + 28N b 1 - 1 ' 2 ir b a r a \ ^d 2 x 

Jr B 

- I [NSa ab D a r b + {2SN + Na- x 5a)D a r a + <j- x 8ar a D a N] ^d 2 x. (66) 

Jr B 

The indices of the perturbed quantities are neither raised nor lowered by a metric. We now assume that a gauge 
may be chosen with r a N a = on the boundary. Under some circumstances this might not be possible as recently 
emphasized by Hayward ]l0| ; but for the applications which we currently envision — namely a static boundary in 
the weak-field zone — this is possible. Thus the first term in the first surface integral of Eq. (|6^) vanishes, and a 
comparison of the coefficients of the variations of N, N a , jah and ir ab with the Einstein equations (^), @, and 
(||), respectively, shows that H\ is a Hamiltonian for the vacuum Einstein equations when the boundary conditions 
are that N, N a and <j ab are given functions of time on the surface at so that the remaining variations within the 
surface terms all vanish. 

A tradition in introductory electricity and magnetism is to consider solutions to Maxwell's equations inside a "box" , 
with the boundary conditions, say, of the potential being given on the surface of the box. The Hamiltonian Hi and 
the boundary conditions of the previous paragraph are appropriate for the gravitational analog of this tradition. 

We have particular interest in periodic geometries, when all of N, N a , -f ab and Tr ab are periodic in t with period 
T = 2t7/Nqloq. It is most convenient to consider the geometrical quantities to be functions of the spatial coordinates 
and a dimensionless time coordinate t = 2-7rf/T; and T is left as unknown and yet to be determined. 

Along with the periodicity, the definition of the time average of a dynamical quantity, say Hi , is 

(Hi) = T- 1 [ Hxdt = ±- [ Hi dr. (67) 
Jo 2vr J Q 

The action used in Maupertuis's principle of least action plays an important role here. This is referred to as the 
M-action A defined by 

167rA= / [ T7 ab Ca a bd 3 xdt 







■2 k 



Tr ab dj ab /dTd 6 xdT. (68) 



2tt 

SNAfyfy + 25N a M a ^ - 5j ab (£ t n ab - V ab ) + 5ir ab (£ aab - Q ah 



d A xdr 



It follows from Eqs. fl66J) and fl6J) that 

+16ttSA/T+— [ I 2SN b -f- 1/2 TT a b r a ^d 2 xdT 
2tt J J rB 

I / [N5a ab D a r b + (2SN + N<j- 1 5a)D a r a + a^Sa^DaN] ^d 2 xdr. (69) 
2tt Jo Jr B 

Eq. ([39j) reveals the use of (H\) in a variational principle. Consider the class of (A, N a ,^ ab , Tr ab ) restricted to periodic 
functions of r, with period 2n, with a specific, predetermined value for A and which satisfy boundary conditions 
at tb specified as particular periodic functions of r for N, N a and <r ab . Then (Hi), considered a functional of 
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the (N, N a , 7 a h, ir ab ) in this class, is an extremum under arbitrary infinitesimal variations of the (N, N a , 7 a fc, ir ab ), 
remaining in the class, if and only if the (N, N a ,j a t,, n ab ) satisfy the Einstein equations. 

Furthermore, the unknown T can be found from the variational principle as well. If only A is changed by a small 
amount, AA, and the variational principle reapplied to find a corresponding change A(Hi), then Eq. (69) shows that 

T = AA/A(Hi). (70) 



IV. THE VARIATIONAL PRINCIPLE EXTENDED TO THE WEAK-FIELD ZONE 

The importance of boundary terms in the Hamiltonian of general relativity was emphasized by Regge and Teitelboim 
jTlJ — different boundary conditions necessitate the inclusion of different boundary integrals in the Hamiltonian. For 
data consisting of N, N a and a ab on the boundary, Hi is the appropriate Hamiltonian. Instead the data could consist 
of, essentially, the derivatives of some of these quantities normal to the boundary, in which case the appropriate 



Hamiltonian would be that of Eq. ( |64J ) but without the boundary integrals at tb- The variational principle of §111 
is closely related to an action principle, which is usually the source of definitions for quantities such as mass and 
angular momentum as particular boundary integrals. Thus, for each different choice of a description of the data at 
the boundary, similar analysis results in a variational principle involving different appropriate boundary integrals, 
different definitions of quantities similar to the mass and angular momentum and a different (H). For the periodic, 
radiative geometries of this paper the amplitudes and phases of the multipole moments (rather than the fields and 
their normal derivatives) are the convenient independent boundary data. And this choice determines the boundary 
integrals of the variational principle and the quantities which correspond to mass and angular momentum. 

In this section we discuss quantities similar to the mass, angular momentum and amplitude and phase of gravita- 
tional radiation, all in the context of periodic solutions of the Einstein equations. Now, these quantities are generally 
not well defined except at spatial or null infinity; and even there technical difficulties persist except under the best of 
circumstances. In addition our analysis is limited to a bounded region of space-time which specifically excludes any 
infinity. It is not surprising then that our approach does not directly involve an elegant, gauge invariant definition of 
any of these quantities. 

Thus we forsake explicit gauge invariance for a precise variational principle involving the multipole moments of 
the geometry. While the physical interpretation of these moments may be suspect, they have the virtues of being 
precisely defined in the context of the chosen gauge and straightforward to calculate. Furthermore, after the fact, a 
periodic geometry can be analyzed and if the amplitude of gravitational radiation (as defined in this gauge dependent 
manner) is sufficiently small, then the geometry near the boundary is of the form of the general linearized solution of 
§ 11 B[ And the mass monopole, X°, and current dipole, S® are reminiscent of mass and angular momentum. 



A. Specifcation of Boundary Data 

The first task of this section is to describe the data on the boundary in terms of the amplitudes , SJL) and 
phases (#/,", Of^) of the multipole moments. The physical metric near the boundary is not necessarily flat and may, in 
fact, be far from flat. None the less, in the vicinity of the boundary we sometimes use a spatial, Cartesian coordinate 
system, (x, y, z) along with the Euclidean tensors hij and S\ defined in Eqs. ( p8| ) and fl57|), in terms of the multipole 
moments; also, the two-sphere at constant r has a unit normal rg when embedded in flat space and r a when embedded 



in the geometry described by the metric ^ a b- The notation of §1IB is employed here except that the symbol r used 
as a tensor index on one of the Cartesian tensors denotes the implied contraction of the index with rg . 

The geometry is still assumed to be periodic, so all tensors are periodic functions of the dimensionless time coordinate 
t with period 2tt. Also the frame in which the geometry is periodic is assumed to be uniformly rotating with respect 
to the frame of reference tied to (x, y, z) with angular velocity Q about the z-axis. 

A prescription for uniquely specifying data on the boundary follows: Choose Nq and tb, once and for all — these 
are never changed. Now choose values of T, SI, I™ m , Of™, S™ m , and 6f™. In terms of these chosen values let each of the 
follo wing geometrical quantities have a value exactly equal to that given by just the linearized theory as described in 
Eqs. (|f to (H): yfd, N 2 ^, ^a ab , N a , N(a c a a*D c r d - \<J ab D c r c ), N~ 1 D a r a , r a D a N and 7 ~ 1/2 ^A&\/^ Also, let 
S a and $ a represent the contravariant tensors in generic coordinates whose components in the Cartesian coordinates 
are given in Eqs. (|57|) and (|5^) respectively. Note that each of these quantities is independent of the others. For 
example, the determinant of y/acr ab is unity and does not depend upon y/a. And N (o^o^D c r & ~ ^a a bD c r c ) is trace 
free and independent of D a r a . These specific choices for data on the boundary are made so that the boundary integrals 
in Eq. ( |73| ) below are precisely equal to what would be expected through second order in an expansion in powers of 
the deviation from flat space. 
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B. The Variational Principle 



The next step toward a variational principle is the definition 



16ttH 2 



N 



*t H TT b ~ TT TT ab 

7 V2 



+ 2N a D b (TT b a / 1 1 / 2 )\^d 3 x 
+ I 2N S a -f- 1/2 ir a b r b V^d 2 x- I 2N<r ab D a r b ^ d 2 x. 

Jr B Jr R 



(71) 



An arbitrary, infinitesimal variation of N, N a , jab and ir ab , which holds fixed the location of the boundary, results in 
an infinitesimal change of (H 2 ) which is similar to the change in (H±) in Eq. ( |69| ) and is simplified by the substitutions 
of A, defined in Eq. (|6£|), and of a quantity, J, which reduces to the current dipole moment for the specific boundary 
data given; 



8ttJ = --L / I <5> a 1 - 1 l 2 n a b r b ^d 2 xdT 
2tt Jo J rn 



= 8nS [ 



(72) 



Then 



l6n5(H 2 ) = — - / / [E Eqns] d 3 x dr + 16tt8A/T + 16^5 J 
2tt Jq J 

S(N a S a ) 1 - 1/2 TT a b r b ^d 2 x dr 



1 

7T 

1 
2^ 



Jr B 
2tt 



Jr B 



5(^a ab )N(D a r b - -a ab D c r c ) + 5(N 2 ^fa)N~ 1 D a r a 



+2S(y/a)r a D a N] d 2 xdr. 



(73) 



The symbol [E Eqns] is an abbreviation for the first integrand in Eq. ( p9[ ) and vanishes when the Einstein equations 
are satisfied. 

The two surface integrals in Eq. (|73[) are examined in the Appendix where it is shown that 



5(N S a ) 1 - 1/2 Ti a b r b ^d 2 xdT 



1 

'2tt 



Jr B 
2n 



Jr B 
2tt 



S(V^a ab )N(D a r b - -a ab D c r c ) + S(N 2 ^/cr)N~ 1 D a r a + 25(^)r a D a N 



d x dr. 



= 5 ^1 r l [2S a 7 -V\ ab r b - 2 2L Da r a WtdPxdT 

1 5 1 

[-h ij r^V k hij - r 0p Vi(h i:i ao jq h qp ) + —h rr h rr - -/ly/iVo + Ir l VJ 



47T 



271 



Jr B 



+S l S l /r-S l riV J S l +^rrgV ' p (fcy)rgv '*) - ^rgV p (rrgV,^*)] ^d 2 x dr 



^(nu -—)S 



Sir 



Jr B 



\2J+1 



KbW) 2t+1 (* + !)(* + 2) fn 32(o; mtl ) 2 '+ 1 Z(/ + !)(/ + 2) Sn 

{l\y^l(l — l) hn lm lm ~> ^ _|_ J^|j2^ ]^ ^Im^lm lm 



(74) 



The different terms on the right hand side of this equation come from expressions in the Appendix. Inside the 
total-<5, the first term is built into expression (All); the next t hree are th e seco nd integral of expression (A14); the 
three terms involving the stationary moments are expressions ( A13| ) and (AlC); and the last two terms come from 
expression (A20). The part of Eq. ( [74] ) proportional to (nu>o — mCl/No) is the second integral of expr ession ( A2C| ). 
And the last summation in Eq. (j74|), which is proportional to the variation of the phases, is expression ( A16). 
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Consideration of Eqs. (f73j) and (74) leads to natural re-definitions of J and A which move the parts of Eq. ( [74] ) 
proportional to 167rf2 into 5J and the parts proportional to 8NqUq(= IQtt/T) into SA. These changes dramatically 
simplify Eq. ([73ft to the result in Eq. (F9J) below. Thus the effective angular momentum, J*, is defined by 



2tt Jo Jt b 



(75) 



And the effective M-action, A*, is defined by 



16tivL 



ir ab d~/ ab /dTd 6 xdT 



E 



2u> r , 



KV#-) r ?V,(C') - h^*r^ p (rr^ q h^)]^d 2 x. 



(76) 



Similarly it seems natural to redefine (H^) by having it absorb the terms inside the total-5 in Eq. (|74|). But, instead, 
we introduce the effective mass which is a boundary integral independent of Nq: 

1 f 27T f 2N 1 f 27r f 

167rm* ee -— / * ——D a r a y/ad 2 xdT+— / ® 2\7 y^Ja^ d 2 x dr 

2n Jo Jt b N 27T 7 J Tb 

+ — f I 2S a -f-^ 2 Tr a b r b ^d 2 xdT 
2vr J J rB 



1 

47T 



2n 



Jr B 



1 5 1 

-h lJ TqV k hij - r p a V iQi 10 <j j q h q p ) + —h rr h rr - -hijh 1 ^ 



+Ir a V l I + S l S l /r-S t r 3 V J S l 
+l<V p (h lJ )r Q >V q (hv) - i/H,-rgV p (rr«V g /i«) 



<tq d x dr. 



(77) 



The variation of the second integral is identically zero; it is included so that m» has the expected value for, say, the 
Schwarzschild geometry. 

Now a variational principle for is based upon 



lQnN a m* = — — 

Z7T 



2tt 



N 



if- 

7 \2 



" + _ TT 71 ^ ^ ~ n n ab 



+ 2N a D b {n b a / 1 1 / 2 ) \ ^d 3 xdr 



+N x [right hand side of Eq. (77)] 
And the variation of this equation yields 



(78) 



167dV <Sm* = ^- J [E Eqns] d 3 x dr + 16ir5A*/T + WnSlSJ* 



Iran 



8(u. 



inn) 2L + 1 (l + + 2) „ r n*r/)/n 

— HmHm ou lm 



(l!) 2 l(l-l) 



32(u mn ) 2l+1 l(l + l)(l + 2) Qn QmxoSn 



[(i + l)!] 2 ('-l) 



D lm°lm 0a lm 



(79) 



To apply the variational principle: choose specific values for J», A* and each of the 9f™ and and consider 
the class of periodic geometries described by (N, N a ,j abl ir ab ) which have these specific values. Eq. ([79]) shows that 
m* , evaluated by Eq. ( f78| ) , is an extremum for a member of this class if and only if the geometry is a solution of the 
Einstein equations. 

As might be anticipated, the coefficients of <S#j"/27r and <56>^"/27r are just 16tt times the energy in one wavelength 
of the gravitational waves in the wave zone as derived from Thome's Eq. (4.16') §. 
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C. The Interpretations of A*, J*, and m» 



In this section we show that when the boundary is in the weak-field zone, each of the effective quantities is 
independent of the exact location of the boundary through terms of second order in the deviation of the geometry 
from fiat space. Thus, for each of the equations of this section, equality is understood to hold only through terms of 
second order. 

The quantity A* is defined as a sum of a volume integral and a surface integral in Eq. ([76|). And the difference in 
A* corresponding to two different boundary surfaces, at r\ and r%, is 

16ttAA, ee 16tt(A 2 * - = / / ir ab dj ab /dT d 3 x dr 

Jo Jri 

-E^-f/-/) KVp^JrgV^^J-^rgVpKV,^")]^^!. (80) 

m,n mn V™ Jr ^ 

Through terms of second order in the weak-field zone this is 

Trmj mn h% n h% n *^d 3 x 

7/1,71 1 

Tin r T2 

-Eo / Vp[rVp(hTKV q (h%n-h% n *VArr%V q h%n]V^d 3 x. (81) 

And by use of the wave equation, Eq. jo^), for h™ n this simplifies to 

i6^a^„ = / 2 i^^^'^rf 3 . - / 2 2^ m ^ n/i r*v^^ 2 ^ 

= 0. (82) 

Thus the effective M-action is independent of the surface over which it is evaluated, through second order in the 
deviation from flat space. 

Similarly, the difference in effective angular momentum evaluated at two different surfaces in the weak-field zone is 
8ttAJ* ee 8tt(J 2 » - Ju) = --^ J (J> - j> ^ <f> b -f- 1/2 ir b a r aV ^d 2 xdT 

+ E^ (<f f ) KVp^KV^^J-^rgVpKV,^)]^^. (83) 

m.n mn \Jr 2 Jri/ 

The first integral may be rewritten as 



7T I I ^ D^-^TT^^fd^dT 
JO Jri 



1 

2^ 



Jri 



' 1 ' Da^- 1 ' 2 ^)^ + 1 - 1 / 2 n ab D a ^ b \ ^fd 3 xdr. (84) 

The first term vanishes from the momentum constraint, Eq. (p), and the second is easily written in terms of the 
deviation from flat space. Also use of the wave equation, Eq. (|63), for h%- n in Eq. (83) results in 

= 0. (85) 

Thus the effective angular momentum is independent of the surface over which it is evaluated in the weak-field zone, 
through second order in the deviation from flat space. At first order J* is just the first term in Eq. (f75|), which 
evaluates to 5° , the current dipole moment as defined by the multipole structure of the geometry in the weak- field 
zone. Thus we have the following interpretation of the effective angular momentum: In the weak-field zone, at linear 
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order, </* is the constant current dipole moment, 5°. But when the second order corrections are added to the linearized 
theory, then S®, has a small contribution from the standing gravitational waves which is linear in r. However, J* 
does not change at second order. Thus, through second order in the deviation from flat space J* is the current dipole 
moment of the source alone without a contribution from the standing gravitational waves. 

In a similar fashion the difference in effective mass evaluated at two different surfaces, Am*, is defined from Eq. (|77|). 

Now, Am* is quadratic in the deviation of the geometry from flat space — the linear contribution vanishes from 
Eq. (|7^) and the fact that A m* is zero for flat space. It is straightforward but an exceedingly tedious task to substitute 
the linearized geometry of § [IB into the definition of Am*. In the later stages the details of the substitution and 
reduction are similar to that for AA* and A J*. And the result is that Am*, too, vanishes through second order in 
the deviation from flat space in the weak-field zone. 

Thus the effective mass is independent of the location of the surface over which it is evaluated in the weak-field 
zone through second order in the deviation of the geometry from flat space. At first order the effective mass is the 
first two terms in Eq. ( |77| ) which evaluates to X°, the mass monopole moment of the linearized geometry of § II B . 
Thus we have the following intrepretation of the effective mass: In the weak-field zone, at linear order, m* is the 
constant mass monopole of the linearized geometry. But, when the second order corrections are included, then the 
mass monopole has a small contribution, linear in r, from the standing gravitational waves. However, m* does not 
change through second order in the deviation from flat space. Thus, to the extent that it is possible to define such a 
quantity, m* is the mass monopole of the source alone without a contribution from the standing gravitational waves. 



V. CONCLUSIONS 



Paper I showed that the periodic solutions of the Einstein equations, with standing waves, give valuable information 
about physically realistic systems, with outgoing radiation. And the variational principle of Eq. ( f78| ) is a valuable 
tool for studying the periodic solutions. In a preliminary application in Paper I, the variational principle has been 
relatively easy to apply in the study of close orbits of equal mass black holes, without regard to the gravitational 
radiation. The next paper in this series will use more sophisticated trial geometries which form a complete set via an 
infinite sequence of possible parameters. In particular these geometries will clearly contain gravitational radiation. 

For specific values of J*, A*, 9^ and Qf£ and an approximate solution to the Einstein equations, accurate to order 
5, the variational principle directly gives an estimate of m* which is accurate to order 5 2 . But the variational principle 
can be used to estimate ft, uiq, llj^] 2 and |<S]^J 2 to order S 2 as well. For example, to find an accurate estimate of ft 
just apply the variational principle a second time with a slightly different choice for the effective angular momentum, 
say a value of J* + A J*. The resulting value of the effective mass will change by an amount, say, Am*; then Eq. ( |79| ) 
shows that 

ft = 7V Am*/AJ* +0(5 2 ) (86) 

accurately estimates ft. The other interesting quantities describing the geometry can be found similarly. 

Perhaps the most interesting results of this paper are the definitions of the effective quantities — Eq. ( f77| ) for the 
effective mass, Eq. ( [75] ) for the effective angular momentum, and Eq. ([76]) for the effective M-action. These definitions 
capture the concept of the mass monopole, current dipole and M-action of the source without contribution from the 
gravitational radiation in the weak-field zone. They are independent of the location of the boundary, through second 
order in the deviations from flat space, as long as the boundary is in the weak-field zone. A similar analysis has been 
made of a toy problem of a mass, confined to oscillations along the z axis under the influence of a generic potential, 
and attached to a semi-infinitely long string stretching out along the x axis. For that case the analogous variational 
principle is for the average energy of the mass alone — the energy of the string does not contribute. 

To add matter to this analysis it is simple to include matter terms in Hi. And to allow for black hole sources 
Eq. (fr8|) for m* still provides the variational principle. But the volume integral has boundaries in the weak-field zone 
surrounding the spatial infinity within each black hole as discussed in Paper I. Then the variations bring in additional 
surface terms, resembling the ones in Eq. J79]), but evaluated at these additional boundaries. These new surface terms 
correspond to the effective mass and angular momentum and the phase of the radiation in the weak-field zone within 
the holes. 
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APPENDIX A: BOUNDARY INTEGRALS IN THE WEAK-FIELD ZONE 



The detailed analysis leading to Eq. ( |74| ) is relegated to this app endix . 

The geometry on the boundary is described at the beginning of §IV A. In particular on the boundary a number of 
quantities are exactly equal to their linearized counterparts. Thus, just on the boundary 



\fo = + I - ~h rr ), 



(Al) 



N 2 V^ = N$^(l - -hrr), 



\foo ah = \foo[o^ h {l — -h rr ) ~ CTQ C hcdCro b ], 



(A2) 
(A3) 



iV a = N S a + fi$ a , 



(A4) 



12 1 1 

2N{a c a a c b l D c r d - -a ab D c r c ) = N [--a L Qa h cd a^ b + -a aab a^ d h cd + a L Qa a^ b r e V e h cd - -a Qab (j c Q d r e V e h cd 

-^Oa^Ob^chde - <7Q a aQ b r e V d h ce + CToabCrgVV c hde], 



(A5) 



N- 1 D a r a = N^ + ^1L + r gv a (7 + \h„)l 



(A6) 



r a D n N = -_^rgV„J 



No. 
2 



and 



l- ll2 *a b r b ^ 



cr dh ar 



2N n dt 



vor V {a S b) . 



(A7) 



(A8) 



Now, consider the term proportional to the stationary SS a on the left hand side of Eq. (|74|). It has no contribution 
from time dependent terms: all time dependent terms are multiplied by stationary terms and drop out when integrated 
over a full period of r. With application of Eqs. (0, [42, 46) and use of the periodicity, the terms proportional to SS a 
in Eq. (^) contribute 

^ f I 25S a j- 1/2 ir a b r by ^d 2 xdT = / I SS^WiSj + W j S i y Qy /^d 2 xdT 

2tt Jo J rB 2tt J J Tb 

f I r- 1 S l S i ^d 2 xdr f I 6S i (V j S i y Qy /^ ~d 2 xdT. (A9) 

JO Jr B J 27T Jo J r 



= ^6 
An 



If there are no stationary sources at infinity, then S* 1 is of the form given in Eq. @ and orthogonality of the y% 
implies that the second term above is symmetric in S l and 5S l . Thus the right hand side of Eq. (A9) is 



N . 
Air 



f <f> (S^i/r - S i riV j S i )^d 2 xdT 

JO Jr B 



(A10) 



The other surface integral on the left hand side of Eq. ( |74| ) is closely related to 

,-2tx 



1 

"to Jo 



5{^a ab )N{D a r b - ^a ab D c r c ) + Si^y/FjN' 1 D a r a + 25{^)r a D a N 



d 2 x dr 



AT p&Tt r i 

2^ Jo J l^^khi^V^^xdr, 



(All) 
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The part of expression (All) depending upon the stationary moments is 



N f 27T f 



(A12) 



And if there are no sources at infinity, then Eq. ( p6| ) and the orthogonality of the 3^™ imply that this is symmetric 
in I and SI and equal to 



No 

47T 



/ f Ir l V 'iIy/ao~d 2 xdT 

JO Jr B 



(A13) 



Expressions (A10) and ( A13| ) are the contributions from the stationary moments to the right hand side of Eq. ( |74|) . 

The radiative terms in expression (All) follow with the substitutions from Eqs. (Al) to ( |A8| ) along with Eqs. (|43|). 
(p5|) and use of Stoke's Theorem; the result is 



jy / f 1 1 

~r~ I & (-8h^r^V k h ll --h^r^\7 k 8h lJ ) y ^d 2 xdT 
4tt J J rB 2 2 

{/" f r ^ r0p ^ <J ° jqhPq ^ + -^,h rr h rr - ^h i] h l k a ik ]^/a^d 2 xdT^ . 



N 0x 
Air 



(A14) 



The focus remains on the first integral of expression ( A14 ); the second integral appears directly in Eq. (|74|). 

The quantity hij is decomposed into its m, n components as in Eq. (|62|), where h mn is an exact solution to the 
fiat space wave equation ([33|) with frequency oj mn and dependence upon I™ m , Of™, Sf m , 0f£ and also upon SI and u a 
through the dependence upon u) mn . Thus 



Shi 



E 



6h% n 

XTn 01 lm 
U1 lm 



5h"" 1 , Sh™ n 5h^ n q 

l 3 xqln I l 3 x an , l J xnSn 

SS n lm ~ ; 



qln 

7 hn 



I in 



56? n 

lm 



>lr, 



SLO r , 



(A15) 



Now, the substitution of Eq. ( |A15| ) into the first integral of expression ( |A14| ) and the use of orthogonality properties 
of the y l £ leads to a sum of terms each involving a single choice of I, m, n and either I™ m or Sf m . The resulting terms 
that are proportional to SI 



and to 66f™ are easy to evaluate because the coeffi cient s of 61% 



and (5(9/™ in Eq. (A15) 



are also solutions to Eq. (J63|); and the terms in the first integral of expression (A14), which involve these, form the 
Wronskian of two solutions to Eq. (|63|). In fact, the term involving 81™ m is just proportional to the Wronskian of two 
dependent solutions of the same linear equation and, hence, vanishes. The surface integral of the Wronskian involving 
X ^Ln i s independent of the surface over which it is evaluated; and it is easiest to evaluate in the wave zone. The terms 



in the first integral of expression ( A14 ) involving 8S\ L m and 80f£ may be handled in a completely analogous manner. 
This results in 



N< 



Imn 



Q*+ 1 (Z + l)(Z + 2) 

(W-i) 



1 lm 1 lm c 



32(u, mn ) 2l+1 l(l + l)(l + 2) 
[(f + l)!P(Z-l) 



qn Qn*rn 



Sn 



(A16) 



for the J-phase contribution to the first integ ral of expression ( A14 ); this appear s dir ectly in Eq. (|74|). 

The term in the first integral of expression (A14) (after substitution from Eq. ( A15| )) which involves 8u) mn is rather 
more complicated. From Eq. (^3|) it is clear that the product r/i™ n may be written so that r and oj mn only appear in 
the combination ruj mni except for a possible overall amplitude dependence upon u> mn . Thus 



5uj n 



8h r - 



1 r 
OUJ r , 



o{ruj mn ) 



and 8(rh mn ) /8(rcj mn ) is closely related to just the r derivative of rh^ n , so that 



(A17) 



8^ 



Sh% n 

0UJ 7nn 



8uJ mn 9 



mn \ 

dr v " ij ' 



Slu 



mn k 



r fc V fe (WC"). 



(A18) 



Now, the 8tu mn contribution to the first integral of expression (A14) is (let 
^w TOn rgV k (rh™l n ) /u>mn, sum over m and n and do the integral over r) 



h% n * and Shij 
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TE/ ^KV p (r/CKV./C* - ^ V p[ r oV 9 (C)]}^^^ (A19) 

„ J t ft LUmn 



lr B 

The term in braces simplifies, and a 5-variation by parts and average over r result in 



'8tt 



/ H ZJ 

'r B 



N ° Jr £ 

-ft™VSV P KV,ft™)]^d 2 I( ir} ; (A20) 
this appears directly in Eq. (|74|). 



'2 , 
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